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We re-interprete the microcanonical conditions in the 
quantum domain as constraints for the interaction of the "gas- 
subsystem" under consideration and its environment ("con- 
tainer"). The time-average of a purity-measure is found to 
equal the average over the respective path in Hilbert-space. 
We then show that for typical (degenerate or non-degenerate) 
thermodynamical systems almost all states within the allowed 
region of Hilbert-space have a local von Neumann-entropy S 
close to the maximum and a purity P close to its minimum, re- 
spectively. Typically thermodynamical systems should there- 
fore obey the second law. 



The second law has been formulated in a number of 
different ways H: According to Clausius' postulate heat 
cannot flow spontaneously from a colder to a hotter sys- 
tem. Thomson's formulation reads: "It is impossible to 
construct a perpetuum mobile of the second kind." Such 
rules define, essentially, some sort of irreversibility, i.e. 
the existence of a state function usually called entropy, 
S, which can only increase in closed systems. 
The second law is arguably one of the most fundamental 
and far-reaching laws of physics; nevertheless its origin 
remains puzzling. 

On the macroscopic level the second law enables us to 
calculate state equations from the requirement that en- 
tropy has to be maximized under the condition of given 
extensive variables in order to get a thermodynamical 
potential as a function of those. 

On the microscopic level irreversibility comes in conflict 
with the notorious reversibility of all fundamental phys- 
ical laws. It has been the challenge for statistical me- 
chanics to reconcile dS/dt > with the underlying mi- 
croscopic dynamics. 

Irreversibility in classical mechanics is conventionally in- 
troduced via two different schemes (cf. ra): The Boltz- 
mann approach, based on the hypothesis of molecu- 
lar chaos (Stosszahlansatz) and the Gibbs ensemble ap- 
proach, based on the hypothesis of quasi-ergodicity. Both 
these attempts have to acknowledge some additional as- 
sumptions which do not follow from the underlying mi- 
croscopic laws. 

This fact has led to continuing controversies as to 
whether this situation can be the last word. Several 
researchers suggested quantum mechanics as a possible 
remedy. 

L. D. Landau and E. M. Lifshitz H expected that quan- 
tum measurements - making a difference between future 



and past - should be responsible for the entropy increase 
according to the second law. This assertion has never 
been proved, though. 

E. Schrodinger || argued that the conventional canoni- 
cal distributions of thermodynamics can be obtained also 
for a system described by a single wavefunction, if one 
invokes "the old crux of molecular disorder" . 
J. von Neumann || was able to show that a non- 
degenerate system should, indeed, obey a quantum- 
mechanical version of ergodicity under fairly weak con- 
ditions. In addition he felt obliged to introduce "macro- 
scopic observers" (coarse graining) in order to come into 
contact with standard thermodynamics. 
This latter aspect has then been taken up also by W. 
Pauli and M. Fierz ||. These authors insisted that the 
2nd law should be explained without reference to exter- 
nal perturbations like those induced by quantum mea- 
surements. 

G. Lindblad j?J observed that the entropy of a multi- 
partite system defined as the sum of the respective par- 
tial entropies of the subsystems would tend to increase 
due to the neglect of correlations (entanglement). This 
fact could be taken as a quantum mechanical justifica- 
tion of Boltzmann's Stosszahlansatz. 
W. Zurek || and his coworkes have discussed in great de- 
tail how the interaction of a quantum system with its en- 
vironment may induce quasi-classical behavior in the for- 
mer ("environment-induced superselection rules"). They 
argued that the second law should eventually result from 
the impossibility of isolating macroscopic systems from 
their environment Q. 

Many of these ideas are pertinent to our present investi- 
gation. It is the purpose of this Letter to show that under 
(appropriately redefined) microcanonical conditions the 
second law of thermodynamics follows from a quantum 
mechanical analysis of the total system partitioned into 
the object and the environment. 

A purity measure that is formulated within standard 
quantum mechanics but has also an interpretation within 
thermodynamics is the von Neumann-entropy 



S(p) := —kTi {p In p] . 
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If this entropy is zero, the system is in a pure state, if it 
takes on its maximum value S max — khi(N), where N is 
the system size (number of accessible states) , the system 
is in the maximally mixed state. Another measure is the 
"purity" P: 



P(p) :-Tr{p 2 } 
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If P takes on its maximum value 1 the system is in a pure 
state, if P takes on its minimum value the system is 
in the maximally mixed state. 

For those extreme cases the two measures uniquely map 
onto each other. For general cases, however, this does not 
hold true anymore. Nevertheless, states with P(p) ~ 1 
will have S(p) « and states with P(p) « P m in will have 
S(p) « S m ax- In the following we mainly consider P and 
get back to S in the end. 

Using the von Neumann-equation for the density opera- 
tor of the total system it can easily be shown that the 
entropy and the purity of a closed system that does not 
interact with any other system are conserved. 
This fact might be considered a contradiction to the sec- 
ond law which demands that entropy should be maxi- 
mized during the evolution of a closed system. 
If the full system is being regarded as divided into two 
subsystems (I and II), the reduced density operators are: 



Pi = Tr /7 {p} , p u = Tr 7 {p} 
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Entropy and purity may be defined for each subsystem 
as before, using pi (pu) instead of p. 
Since those are not the density operators that appear 
in the von Neumann-equation, Si (Sn) and Pi (Pu) de- 
fined on basis of these operators are no longer conserved. 
Those are the quantities we are going to examine. 
Although the following ideas apply to all sorts of subsys- 
tems, we want to refer to the system of which the entropy 
is to be calculated as the "gas-system" , g, and all the sur- 
rounding as the "container-system" , c. 
The full Hamiltonian is now divided according to the 
same scheme: 



H =: L a + L c + W 
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L g describes the energies arising from the gas parti- 
cles alone, including their mutual interactions. L c de- 
scribes the corresponding energies of the container par- 
ticles alone. W describes the interaction terms that de- 
pend on both, the coordinates of gas particles and the 
container particles. These are here basically the "wall" 
interactions that keep the gas particles inside the con- 
tainer. 

The energy eigenstates of a free gas are unbound and con- 
tinuous. It is far more convenient to have bound states 
for the "separate" systems. Thus we modify the Hamil- 
tonian formally in the following way: 



H=:L' g + L c 



W' 



with 



L' :=L g + V 



and W' := W — V 
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where V models the mean effect of the container on the 
gas particles. V is an effective potential that only de- 
pends on the coordinates of the gas particles and is chosen 
to minimize W' . Usually the container is simply modelled 



by some "box" potential V, neglecting W' altogether. 
But however small, starting from first principles W' will 
always be present, and represents a coupling. 
We consider the gas-system to be closed in the thermo- 
dynamical sense, i. e., controlled by microcanonical con- 
ditions (E, V,N — const.); this is clearly an idealization 
but can routinely be realized in an approximate way. 
A system that is closed on the macroscopic level (ther- 
modynamically closed) does not need to be closed on the 
microscopic (quantum-) level (i. e. not interacting with 
any other system). 

The fact that no extensive quantities are to be exchanged, 
however, puts constraints on its Hamiltonian, especially 
on the interactions a system can have with its surround- 
ing, in order for it to be thermodynamically closed. 
The strict conservation of particles N is taken into ac- 
count simply by the way the system is partitioned. To 
which accuracy the volume V stays fixed is set by L c . 
Microcanonical conditions then correspond to a box with 
very high, in the limit of the volume V being exactly con- 
served, infinitely high potential walls. 
The condition that no energy is to be exchanged further 
constrains the Hamiltonian. The energy contained in the 
gas is given by: 



J g ■- \-gl 

If this is to be conserved, it follows that 







L' g ,W 



0. 
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Except for these constraints we need not specify W' in 
more detail. 

Based on these commutator relations we find that for any 
energy eigenspace A, B 



Ei<W = Eiv4>)i 2 



(0.9) 



is a conserved quantity, set by the initial state, where 
ip£j B denotes the amplitudes of the degenerate product 

energy eigenstates of L' g + L c ("i" denoting the gas, 
"j" the container part of the product) that are associ- 
ated with the energy eigenvalues E 9 A (Eg) in the gas- 
container-) system. 

Since we want to consider cases here that have zero local 
entropy in the beginning (product states), we get 



E 



i^ B (°)i 2 



>)IW(0)| 2 



p9 pc 
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where P 9 A (P B ) are the probabilities of finding the gas- 
container-) sytem somewhere in the possibly highly de- 
generate subspace characterized by the energy eigenval- 
ues E A (E B ). If no energy is to be exchanged, clearly 
these probabilities have to remain conserved. This is the 
constraint that microcanonical conditions impose on the 
accessible region of Hilbert-space. 
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Although we are interested in P(t — > oo), we start by 
considering the time average of the purity P for reasons 
that will become clear later. 
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P(Mt)))dt 
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Choosing a special parametrization for \ip), we can con- 
vert the time integral into an integral over the trajectory 
generated by the total system's dynamics for given ini- 
tial conditions. Parametrizing \ip) in terms of the real 
and imaginary parts of its amplitudes 



we can write instead of ( J.llj) 



P = 



rlV-(T)) _j_ 



(0.12) 
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where \d\ip)\ denotes the "length" of an infinitesimal step 
along the trajectory in Hilbert-space. 
The advantage of this special parametrization derives 
from the fact that the effective velocity 



u eff 



E 



= ±(m\H 2 \m) (o.i4) 



is constant on each trajectory and thus independent of 
the time t or the special point on the trajectory. Hence, 
the integral fl0.13| ) simplifies to 



P = 



L 



\HT)) 



V>(0)> 



p{{^'ij}) 
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where L is the length of the path. So, the time average 
of P equals the path average along a special trajectory 
in this parametrization of Hilbert-space. 
We are not able to compute this integral for we do not 
know W' in detail, and even if we did, we could never 
hope to solve the Schrodinger-equation for a system with 
about 10 23 de grees of freedom. 

All we want to prove here, is that for typical trajectories 
staying within the region allowed by the microcanonical 
conditions, P is extremely close to its minimum value for 
almost all points within this region. 
We proceed as follows: 

First we calculate P m in which is the smallest possible 
value of P within the allowed region. Then we compute 
the average of P over the total allowed region. If this 
average is close to P m in > we can conclude that P w P m in 
for almost all points within this region, which means for 
almost all P(t), since any distribution with a mean value 
close t o a boundary has to be sharply peaked. 
From ([l^) it follows that the P 9 A remain conserved. Now 
the lowest purity P of any state consistent with this con- 
diton is: 



p V 

L min — / 



9\2 



N% 



(0.16) 



where N A is the degree of degeneracy of E a A . 
To calculate the Hilbert-space average of P denoted as 
< P > we need a para metr ization for ipij , ip' ^ confined 



to the allowed region (3.9) that essentially consists of 
hyperspheres in different parts of the Hilbert-space of 
the total system. The Hilbert-space average can then be 
written as 



< P >= 



/P({^({0 n }),^({0„})})det^n„' 



/det^IL* 



(0.17) 



where <p n is the respective set of parameters and J- is the 
corresponding functional matrix. 

This integral can actually be solved analytically. The 
techniques are essentially the same as those used to cal- 
culate surface areas of hyperspheres in the classical sta- 
tistical analysis of the ideal gas. Since this calculation is 
rather elaborate we do not want to present it in detail 
here, but give and discuss the result: 



< P >= 
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(0.18) 
(P 9 a) 2 ) 
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N A N% + 1 
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Here Ni is the degree of degeneracy of the energy eigen- 



value Eg. 



If the degeneracy of the occupied energy levels is large 
enough so that 
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N%Nl 



1 



N%N C B 
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which shou ld hold true for typical thermodynamical sys- 
tems, (IIS) reduces to 



< P >c 
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(P 9 aY 
N A 
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(P C B? 
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The fi rst sum in this expression is obviously exactly P m in 
( |0.16|) , so that for systems and initial conditions in which 
the second sum is small the allowed region almost only 
consists of states for which P as P m in- The second sum 
will be small if the container system occupies highly de- 
generate states, typical for thermodynamical systems. 
To illustrate this result, we have plotted the relative fre- 
quency of P (see Fig 1.), which we calculated using a 
formula by Lloyd, Pagels [jl0| and Page |ll|] , that applies 
to completely degenerate subsystems only. In this case 
we find from (0.18) 



< P >= 



N° + N C 



N3N C 



1 



(0.21) 



(For this special case the average has also been calculated 
by Lubkin fll.) 
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FIG. 1. Probability density as a function of P = Tr{pf} 
for various nixn2 systems. With increasing size of subsystem 
II the density becomes peaked near the value 

Obviously our average < P > is in perfect agreement 
with the more detailed distribution of P. 
So far we have only shown that, except for negligible 
parts, the Hilbert-space-section in which the trajecto- 
ries can be found under microcanonical conditions, has 
P w P m i n , provided the surrounding system is much 
larger than the considered system. 

To examine under what conditions trajectories will even 
"fill" the whole space they can possibly live in, we con- 
sider the special case of no degeneracy in either subsys- 
tem. Here the evolution of P can be calculated exactly: 



P(t) 



A,B,C,D 



o(Tf[( E AB-EcB+EcD-EAD)t) p9 pC pCJ p 



A r B r C r D 



(0.22) 



where Eu are the energy eigenvalues of the respective 
energy eigenstates. Note that without any interaction 
Ei j = Ei + Ej and P(t) = c onst ., as expected. 
Assuming that all terms in (0.22) are oscillating except 
the ones with A — C or B = D, we get 

A B \ A J \ B ) 

(0.23) 



which is exactly the same result as ( 3. 18 ) with Ni 



1. So, under this assumption the Hilbert-space 



average is exactly equal to the time average. 
In general, however, this specific ergodicity is not needed: 
It suffices that typical quantum trajectories, even though 
starting with P(0) = 1, venture out into the vast Hilbert- 
space regions characterized by P = P m in- Of course, one 
cannot exclude that in special situations there might be 
trajectories that never leave the very tiny region with 
P ~ 1, but these situations become extremely rare as 
the surrounding gets big. 

Finally, we return to the local entropy S. Trying to com- 
pute < S > rather than < P > we get, after some lengthy 
but straightforward perturbative calculations 



< S >» S„ 



{{P%N%})-K 



E 



{P C B? 

Ng 



(0.24) 



where K is a positive function that s cales linearly with 
the system size of the gas system. ( 0.24 ) is valid for 
situations with 
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which is the thcrm odynamical regime in which the sec- 
ond term in (0.24) will be small. (Again, for the special 
case of both subsystems being completely degenerate, our 
results are in perfect agreement with a result by S. Sen 

In conclusion we have shown that the local von Neumann- 
entropy of a considered system will be maximized dur- 
ing its evolution, even if the system is thermodynami- 
cally closed, provided the energy eigenspaces occupied 
by the surrounding are much bigger than the energy 
eigenspaces occupied by the considered system. This is 
typical for thermodynamical systems. We did not need 
the additional assumptions underlying classical deriva- 
tions. Since we considered microcanonical conditions we 
get the maximum entropy as an explicit function of the 
initial energy distribution. This allows for a connection 
with standard thermodynamics, entropy being a thermo- 
dynamical potential. 

Generalizations to canonical conditions and the treat- 
ment of quantum computer systems as specific open sys- 
tems are under way. 
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